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Abstract 

We investigate the noncommutative effect on the non-Gaussianities of pri- 
mordial cosmological perturbation. In the lowest order of string length and 
slow-roll parameter, we find that in the models with small speed of sound the 
noncommutative modifications could be observable if assuming a relatively low 
string scale. In particular, the dominant modification of non-Gaussianity esti- 
mator /nl could reach 0(1) in DBI inflation and K-inflation. The corrections 
are sensitive to the speed of sound and the choice of string length scale. More- 
over the shapes of the corrected non-Gaussianities are distinct from that of 
ordinary ones. 

1 Introduction 

Inflation [Ij is a very successful paradigm of the very early universe. It can naturally 
solve several very tough cosmological problems without fine tuning. Furthermore, it 
predicts a nearly scale invariant Gaussian CMB spectrum, which has been confirmed 
in the experiments [2]. However, one of the problems with inflation is that there are 
too many inflationary models. It is necessary to find the signatures which could 
distinct various models. With the development of precise cosmology, we expect new 
experiments to constrain a large amount of inflation models and to make the paradigm 



more clear. The scalar spectral index and its running, the gravitational wave and non- 
Gaussian component of the primordial fluctuations are among the important probes 
to detect different inflation models. In this paper, we mainly discuss non-Gaussianity. 

In recent years, the non-Gaussianities of primordial perturbation opens another 
window other than power spectrum to study different scenarios of inflation models. 
After systematic analysis [H HJ [5l [6], it is found that non-Gaussianity estimator f^L is 
of order 0(e) in the ordinary "slow-roll" inflation and thus is unobservable. In some 
particular models with non-trivial dynamics, such as ghost inflation [10] and those with 
sufficiently small speed of sound, including DBI model pniT2] and K-inflation [20l [21] . 
the unsuppressed non-Gaussianity is potentially detectable in future experiments [3 
[12] . Furthermore, since the shape of non-Gaussianity is more multiple than that of 
power spectrum, it could be used to distinguish different inflation models through 
the classification of those configuration of k modes that determine the maximum of 
three-point function[8]. 

Non-Gaussianity could also be used to study various trans-Planckian physics pro- 
posals. Deviation from the standard Bunch-Davies vacuum was considered in [TJ [T2] 
where the correction is found to be O(^) where A is the energy scale on which the 
modes are generated. In [9], higher dimensional operators were introduced in the 
infiaton Lagrangian. This modification could enhance the non-Gaussian effect if the 
energy cutoff is not too high; but it is difficult to exceed /atl ~ 1. In this paper, we 
consider another trans-Planckian scenario which is based on noncommutative geom- 
etry. 

In string theory, a promising candidate of quantum gravity, one may use perturba- 
tive string or nonperturbative object D-brane to probe the spacetime geometry. Due 
to the extensive nature of the string, or stringy effect, or strong string interaction, 
the picture of spacetime geometry in string theory could be very different from the 
usual one. Especially, very near the cosmological singularity, the usual concept of 
commutative geometry may break down completely. A better description could be 
noncommutative geometry |22]. in terms of the algebra generated by noncommutative 
coordinates 



One natural way to get noncommutative geometry is to consider the D-brane in the 
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presence of constant NS-NS magnetic B field[23]. In this case, the spatial coordinates 
are noncommutative 

[x\x^] = ie'^ , (1.2) 

where 6^^ depends on the background flux, while space-time is commutative. Similarly 
one can obtain space-time noncommutativity by placing D brane in constant electric 
field, but the theory is no more unitary in this case [211 [251 [SSI [271 128] . 

Another way to realize the noncommutative relation among coordinates, especially 
between space and time, is to start from the stringy uncertainty relation 

AxpAt > ll (1.3) 

Equivalently one may assume [t,Xp\ = ill- This kind of noncommutative relation 
have been applied to the study of inflationary universe |30[ [3T1 [32] . In [30], Bran- 
denberger and Ho started from (II. 3p and discussed the cosmological implications of 
such relation. To keep the background isometry intact, they actually considered the 
noncommutative relation between radial coordinate r and time t, however, as we 
mentioned before, the space-time noncommutative relation may violate the unitarity 
of the theory. 

In this paper, in order to keep unitarity, we choose to set the space-time com- 
ponents 6'°* zero. The cosmological imprint based on this kind of noncommutativity 
has been studied by many authors[3ll [35l [361 EH [38]. Meanwhile, without losing 
generality, we choose a particular frame in which the only non-vanishing space-space 
component of 6'^'^ in comoving coordinate is 

We find the calculation is much more involved if one chooses to keep all possible 
spatial noncommutativity but this won't change the main result we obtained from 
this particular frame. We also require that noncommutativity following from (11.41) 
only dominates at small length scale. A simple form of 9^"^ realizes this condition 
is|38] 

e'' = % (1.5) 

where Is is the string length. In this case, the noncommutativity in physical coordinate 
remains to be a constant, [Xp,^^] = ill, noncommutativity of comoving 
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coordinate is diluted by the expansion of universe. However, with this choice, the 
isotropy of the spacetime is now broken. Nevertheless, if the noncommutativity is 
very small, the breaking of isotropy could be ignorable but the physical implication 
could be observable. In this paper, we will focus on this case. 

In most of the inflation models, the non-Gaussianity is quite small and hard to 
be detected. The noncommutative corrections of non-Gaussianity is even smaller. 
This drives us to work on the models with big non-Gaussianity. We will discuss the 
noncommutative effect in string-inspired DBI inflation model and k-infiation model. 
We find that in both models the noncommutative modification of non-Gaussianity 
estimator can reach 0{1) or even larger with a small speed of sound and a relatively 
low string scale. We also determine the shape of corrected non-Gaussianity, which is 
different from the usual ones. Though the non-Gaussianity estimator could be over- 
shadowed, the shape is a distinctive signature to be probed by the future experiments. 
On the contrary, future experiments could set bound of the noncommutativity and 
be used for crosschecking along with other experiments, like atomic experiments. 

The paper is organized as follows. In section 2, we give a brief review of the ADM 
formalism to discuss the perturbations in a general inflation lagrangian. In section 3, 
we study the noncommutative modifications to general infiation models. In section 
4, we focus on two infiation models with observable non-Gaussianity. In section 5, we 
end with conclusion and discussions. 

2 Perturbations in general inflation models 

Let us start with a general Lagrangian P(X, </)) which can be used to describe a broad 
class of infiation models. The action is of the form as follows. 



where (f) is the infiaton field and X = — i(?^'^c},t09i^0. We set Planck mass Mpi — 
(8nG)^^ = 1 and adopt the metric signature (—1, 1, 1, 1). 

In order to proceed, it is convenient to work in the ADM metric formalism. 





N^dH + hij{dx' + N'dt){dx^ + N^t), 
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where hij = a^(l + 2QSij, with a the scale factor which grows quasi-exponentially 
during inflation and ( the scalar perturbation of metric. In this paper, we do not con- 
sider the tensor perturbation. Since in the ADM formalism and A^* are Lagrangian 
multipliers, the metric in terms of ( can be determined by solving the constraint 
equations of A^ and A^* in a given gauge, rather by solving the Einstein equation, 
which makes the calculation relatively simpler. 

There are two main gauges in which we can calculate the perturbation action. One 
gauge is the comoving gauge, that is 6(f) = 0. In this gauge, ( is the curvature scalar 
on a comoving hypersurface and it directly seeds on the later generation of large scale 
structure and anisotropy of microwave background radiation. It is straightforward to 
find out that ( remains constant after horizon exit. Although the physical meaning 
is manifest within this gauge, it is rather complicated to analyze the order of the 
perturbation action with respect to the slow-roll parameters, which is important in 
determining the magnitude of correlation functions. In this gauge, the third order 
action of perturbation is apparently of order 0(e°), where e represents the slow-roll 
parameter; however after doing a lot of integration by part the action is actually 
of order 0(e^)[3l [7]. The other gauge with ^ = is called uniform density gauge. 
The gauge transformation linking the two gauge is C = ^Scj), where H is the Hubble 
parameter H = ^. One could easily recover the exact order O(e^) of the action by 
doing a gauge transformation of the perturbation action of 50 |3] . In order to calculate 
the correlation function of metric perturbation in this gauge, we first calculate the 
correlation function of 6(f) and then transform it to ( just after horizon exit, which 
remains constant ever since. 

Below we adopt the uniform density gauge ^ = to carry out the calculation. 
Substituting the ADM metric into the action fl2.1l) . we get 

S=^J dtd^xVhN{R^^^ + 2P) + i J dtd'^xVhN^^EijE'^ - E^), (2.3) 

where h = det(/ijj), Eij = \{hij — ViNj — VjNi), E = E^. R^^'^ is the Ricci scalar cal- 
culated in the three-dimensional hypersurface with metric h^j , and V, is the covariant 
differential coefficient defined on the hypersurface. Since A^ and A^* are Lagrangian 
multipliers, we can obtain two constraint equations from them, which are 

7^(3) + 2P - 2A^"2Px(0^ - 2N'4)di(f) + N'N'd^(f)dj(f)) - N-^{EijE'^ - E^) = 0, (2.4) 
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- iVPx(2iV-2090 - j^d4dj4>) + 2{V,{N~^El) - V,{N-^E)) = 0. (2.5) 

We divide (f) into the isotropic background 0o(^) aiid the fluctuation = ^q + V?. 
In order to evaluate the third order action of perturbation, we only need to solve the 
equations of and A^* to the first order of since the second and third order 
of solutions, when substituted into the action, will be multiplied with the first and 
zeroth order Hamiltonian constraint of the action respectively and thus vanish. In 
fact, to calculate the n-th order action of perturbation, one only needs the solutions 
of and A^* to {n — 2)-th order [?]. Following [3], we decompose A^* into two parts 
Ni = Ni + diip, where 9jiV* = and Ni is lowered by hij through A^*. Then we expand 
them to the first order of 

N = l + au Ni = Nl'\ ^ = ^^1, (2.6) 

where ai,N-'^^ and ifji are of order 0{(p). Substituting (12.61) into (12. 4p and (12.51) . and 
solving them to 0{{p), we obtain the solution, 

«i = ^V^, Ar« = 0, (2.7) 

-mP^xM^ - (2Pxo0o + 2P^x,xM) , (2.8) 

where the subindex represents the background value. 

It will be more succinct to express the solutions using the slow-roll parameter e 
and the "speed of sound" c^, 

i^^x^P^ 

' dE Pxo + 2XoP 
where E = 2XP x — P is the energy of infiaton field. Keeping with the lowest order 

of slow-roll parameter, equations (12. 7p (12. Sp can be written as 

He 

«i = (2.11) 

5Vi = (^ - ^e)v. - (2.12) 
2H 00 
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We omit the subindex 1 in these expressions in the following section for simplicity. 

Substituting fl2.11l) and (12.121) into (12.31) and expanding it to the second order of 
perturbation, we attain the free field action of fluctuation 

S2 = IJ d'xa^^^' - P,x,{dv?). (2.13) 
The equation of motion is 

// 

(ay;)" + (c^A;2 - — )ay? = 0, (2.14) 
a 

where the prime denotes derivative with respect to conformal time and we have as- 
sumed Pxo and Cs to be time independent. The classical solution is 

y, = -^^==(l + ^c.A:r/)e--^ (2.15) 

where we choose the standard Bunch-Davies vacuum to fix the coefficient. With the 
same procedure one could obtain the third order action of perturbation which has 
been studied thoroughly in \[\ in the comoving gauge. 



3 Noncommutative modification 

In the noncommutative spacetime, the functions are better described by the operators 
in Hilbert space. This is very similar to the case in quantum mechanics, where one 
has a noncommutative phase space. The product of functions may be taken as the 
multiplication of the operators. But an efficient way to define product is by so-called 
Moyal product, whose expansion in curved spacetime gives [38j 

fc=o ■ ^ ^ 

where Z)^ is the covariant derivative coefficient in the curved spacetime. 

The star product could be extended to multiple function situation, at quadratic 
order in 6^'^ , 

/i*---^/n = (l + V'^J]/^^!?^- V^"'^ J2 D;DlD;Di)f,---U (3.2) 

a<b a<b,c<d 

To incorporate the noncommutative effect, we replace the ordinary product in 
the infiaton Lagrangian P{X, 0) with the star product. Suppose that P{X, 0) relies 
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on (f) through function V{(l)) 0, which could be the inflaton potential or the warping 
factor in DBI model. Without losing generality, we simply consider V^(0) of the form 
V{(j)) = 0", n > 1, since one could always do Taylor's expansion of a generic function. 
When apply the star product in the generic Lagrangian, noncommutativity changes 
the forms of dynamic term X = —^d^(j)d^(f) and scalar function as follows 

Z o 

-DiD2D^<PD2DiD,</)) + 0(6^), (3.3) 

and 



V = (h-k ■ ■ ■ -k i 



n(n — 11 



i12n2 jn-2 



(^1^10^2^20 - DiD2(pD2D 



n{n l)(n ^(^i2)2^n-3p^^^^^^^^^^ + ^^^^^^^^^^^ 



24 ' ' ' 

-DiD2(f)Di^D24> - D2Di4>D2^Di4>) 

+0{e^). (3.4) 

We find that the lowest order noncommutative modification term is of order 0{d^). 
We only consider the correction of the lowest order 6"^ in this paper, and denote them 
as 6gX and 6oV respectively, 

6eX = \{e^^fg^^''{D,D^D^<pD2D2D,<P-D^D2D^<pD2D,D,<P), (3.5) 



Tiifi — 1 I 

5eV = ^-—^{e^^fr-\D^D^(l)D2D2<p-D^D2ct>D2D^<P) 

o 

n{n - l){n -^(^12)2^n-3(^^^^^^^^^^^ + ^2^20^10/^1^ 



24 V / - 

-DiD2(l)Di(f)D2(p-D2D^(j)D2(pD^(j)). (3.6) 

One should keep in mind that all the covariant derivatives are calculated within 
the ADM formalism, which makes the evaluation much more involved. In some other 
papers in calculating large unsuppressed non-Gaussianity [lOl [HI [9], because of the 



^If there are more than one such functions, one just needs to include all the corrections from each 
function in p.l4p . 
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particularity of the actions (they could themselves, by doing serial expansion, generate 
cubic terms of perturbation without involving metric corrections, i.e. a, the 
authors ignore the correlation between metric correction and inflaton perturbation 
in the uniform density gauge, which is subleading in slow-roll parameter, and do the 
calculation with isotropic FRW metric. Our calculation turns out to be the same 
situation when considering the noncommutativity correction. 
The change of the inflation action is 



To obtain the third order action of perturbation, one needs to serially expand Px 
and Py around the background value and multiply with the corresponding terms in 
6gX and 6gV which generate overall cubic terms of perturbation. 

To simplify the calculation, we need to pick out the terms of leading order of 
slow-roll parameter. We find in this case, the leading order is 0(e) for ( after gauge 
transformation. We decompose SqX into terms of different order of perturbation. 



where (^0-^)3 represents the cubic terms, etc., and we do not consider the higher order 
of perturbation since we are only going to calculate the three-point function. The 
terms of lowest order of slow-roll parameter in [dgX)^ are those composed of product 
of two inflaton perturbation </? and one metric correction, i.e. a and ip. However, 
we do not need to take them into account to obtain the final third order action of 
perturbation by the reason that they will generate terms of second order in slow- 
roll parameter after doing the gauge transformation to C, which turns out to be the 
subleading terms in our calculation. As for {SeX)^ which is composed of terms of 00, 
it also results in subleading order terms when multiplied with the serially expanding 
terms of Px- In short, we only need {6gX)2 and {6gX)i. 

The terms of lowest order of slow-roll parameter in {d0X)2 are the products of two 




(3.7) 



SgX = {deX)3 + {50X)2 + {6gX)i + {5eX)o, 



(3.8) 
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inflaton perturbation, which are summarized as follows, 

-Hdldi^pdiip + H^dldiipdiip + H^d^diipdiip + d^d^ipd^ip 
-\-a^d2pd2p + a^dipdip — 2Ha^diipdiip — 2Hd'^d20d2(p 
-2Hd^diPdi^ + {Hdfdi^di^ + {Hdf 82^82^ 
+{Hd)'^di(pdi(p — d\(pd2(p + 9i92<^9i92<^ — {HadYpP' 
—Hadpd\p + Hadpd^p + ^Hd^pd^i^ + 2Hd^ipd2p 
+2H'^ad(pdfip + 2H'^adpdl^ - AHdid2pdid2^ 
—AH'^dl^pdl^ + AH"^ 8182^8182'^ + add\ip)(p + add2^p(p 
—a^d^(p(p — (i'(p{2d\i^ + 2d\i^ — adp)), (3.9) 

where i should be summed from 1 to 3. And {5eX)i is 

{5eX)i = --{e^^fHad^oidlp + 8^^ - 2H8lip - 2H8l^p + 2Had(p). (3.10) 
8 

Following the same procedure, we decompose V as 

5eV = (5eV)^ + (5eV)2 + (Wi + {5eV)o. (3.11) 

According to the same reason as that of the case of SqX, we can just consider the 
leading terms in {5gV)2 and {50V)i, which are 

o 

adipdlif — ad(pd2^) , (3-12) 
{SeV)i = - ^(^ ~^) (gi2)20g-2(2(ad)20o(^ - ad0o92(^ - ad0o922</^). (3.13) 

o 

The whole change of third order action of perturbation due to noncommutative 
geometry in leading order of slow-roll parameter can be written as 

SeSs = j d''xVh{P,XoXo{SgX)i{6eX)2 + P,XoXo{SgX)2i6eX), 

1 



+P,Xo^MW2 + -P,XoXoX,{S,X)i{SeX) 



2 



+P,v,x,{S9XUSgV)2 + PyoXo{^X)2{SeV)i), (3.14) 
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where 5gX = X — Xq, and {5gX)i = 000, {SgX)2 = ^(0^ — (dipY), where we have 
picked out the terms with least (pQ to reduce the order of slow-roll parameter. Al- 
though we have picked out the leading order terms, the result is still lengthy. We 
write the result of (13.14^ in Appendix A. To write the result, we define two new 
parameters, 

4 = p , (3.15) 

a = XqXqXqP^XoXqXo- (3.16) 

As we pointed out above, all the terms except those with a coefficient are of order 
0(e) in ( after a gauge transformation (p = These two parameters vanish in 
the general "slow-roll" inflation models, but could be non-trivial in some particular 
models. We will see in some case, they determine the dominance of the correction of 
non-Gaussianity. 

There are some subtle differences between calculating the modification of three- 
point function and two-point function due to noncommutativity. In calculating two- 
point function, one has to solve the equation of motion of perturbation which is in 
general hard to solve with the presence of noncommutative coordinates (for a solv- 
able example, see [M])- In Ref . [30] . the author developed another way to encode 
the noncommutative effect into the power spectrum without solving the equation of 
motion. 

However, even with a solvable equation, we do not need the modified classical 
solution to calculate noncommutative correction of three-point function. The three- 
point function is calculated through 

{C%t)) = -t fdt'{[C'{t),HUt')]) (3.17) 

J to 

in tree level. So, to evaluate the modification of three-point function, which is denoted 
by {C'^{t))0 below, in the lowest order of 6, we divide the Poisson bracket into two 
groups: 

{C'{t))g = -I f dt'{[Cl{t),5gHUt')])-^ fdt'{[3C{t)6gC{t),HUt')]), (3.18) 

where Cc is the commutative solution. Since the primordial Hamiltonian without 
noncommutative modification is of order in slow-roll parameter, thus the second 
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group contributes terms subleading in slow-roll parameter, as we emphasize the lead- 
ing order of the modification part of Hamiltonian is 0(e). While the order of 9 are 9"^ 
in both case, so the leading modification is obtained from the first group of Poisson 
bracket which is summarized in Appendix B. 

Another point deserved mention is that the constraint equations of and A^j, and 
thus the solutions (12.111) (I2.12p obtain corrections of order 9^ in noncommutative coor- 
dinate. However, the same as the above analysis, they contribute terms of subleading 
order of 6* or e either in 5gS oi S and thus we do not need to consider them. 

4 Model testing 

In this section, we evaluate the effect of noncommutativity on the non-Gaussianities 
of perturbation in some particular models. Note all the background field and Hubble 
scale are estimated at horizon crossing, namely at the time about 60 e-foldings before 
ending of inflation. 

In the ordinary "slow-roll" inflation with Lagrangian P = X — V , since 

= 1, cy = 1, a = 0, (4.1) 

it turns out that the modification is zero on the level of first order of slow-roll param- 
eter. Let us turn to the models with significant non-Gaussianity. 

4.1 DBI model 

DBI infiation [miElIISlITl] is motivated by brane infiation scenario flHl [T71 [T8l [T9] 
in warped compactifications. The effective Lagrangian is 

P(X,0) = -/(0)-Vl-2/(0)X + /(0)-^ - V^(0), (4.2) 

where / is the warp factor / = and A depends on fiux number. The value of 
speed of sound c<j, as well as other two parameters c/ (/ substituting for V in (13.151) 
represents warp factor /(</>)) and a are as follows, 

c, = v/l- 02/(0), (4.3) 

Cf = V^l- 02/(0), (4.4) 

c^ = -^0V^(l-0V)-i (4.5) 
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We find c/ coincides with Cg in DBI model, and thus it becomes as small as the speed 
of sound. The leading correction terms with respect to Cg comes from (IB.ip are 



-klk\klkl){k, . h)kl{2A^^ + 120^1^ + 720^) 

- klk'r,kl){k, . h){8^^ 

fci ko + kok-i + fcqfci ki kok'i , 

+24 y + 120^) 

-Q^H^\fimr+ik1r)kM2^ 

-^?fc,^fc,U,^)fc3^(2i, + 6^ + 24^)) 

-^^^mr + (k\r)i^^kikiki^ 

pi 

-2klkl{2^^+Q^))+perm.i (4.6) 

where a and h denote the first and second component of k vector respectively, K = 
ki + k2 + ^3 and perm, denotes all the other terms obtained by rotating the index 

J 4 

(1, 2, 3). Here we have used the relation 0q ~ ^ due to small speed of sound. In this 
phase, we have power spectrum of (p perturbation 

where A'^e ~ 60 is the e-folding between horizon crossing and the end of inflation. 
According to COBE normalization, ^ 23 x 10"^°, we find A ^ 10^^ and 4>o ~ 
10~'^Mpi. Using the limit of small speed of sound, 0o ~ 0q^^A^/^ ~ Mpi. As a matter 
of result, only the first term in (14.61) with the following shape dominates, 

A = -^((^?)'(A:^)'-A:?fc?A;»A;,^)(fci -^2)^:3^(24-1^ 

ki + ko _ kiko. , . . 

+ 120 + 7204;f ) + perm., (4.8) 

We parameterize {C,kXk2Ckz)e as, 

{CkAMe = {27ry6{h + h + fc3)(Pc')' J^-^i' (4-9) 
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Figure 1: The shape of Ai/kik2k^. 



where is the primordial power spectrum [20 



C 

Pt = . (4.10) 

The shape of Ai/kik2k3 as function of Xi = ki/k^ and X2 = k2/k^ are drawn in Fig. [1] 
In drawing the figures, we omit the coefficient of Ai. And we have chosen a particular 
frame in which the three k modes are in x — y plane and k^ is along the x-axis. We 
find this shape is distinct from the shape of ^c[Z|- 

The non-Gaussianity of CMB in the WMAP observations is analyzed by assuming 
the following ansatz 

C = C,-\fNLCl (4.11) 

where C,g represents the Gaussian part of C, and fNL is an estimator of non-Gaussianity. 
The three-point function of C can be factorized as 

= (2vr)^5^(fci + k2 + ^3)(-^/^l(Pc')')^- (4-12) 

Despite the difference between the shape of (14. 6 p and that of f l4.12p . we set ki = 
k2 = ks = kto calculate /ji^^ which represents the size of correction of non-Gaussianity 
deriving from Ai. We have 

fhL = 0.02^, (4.13) 
where the results are evaluated in the particular frame we used above. 
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We find tlie noncommutative correction of non-Gaussianity could be large if Cs <C 
1. For example, for Hubble constant H ~ 10~^Mpi, the string length scale Is ~ 
lO^Mp^^ and the speed of sound Cg ~ 0.1, then flf^ ~ 2. Comparing to the dominant 
Jnl of commutative case((5.10) in [7]), 

/^^^0.32c;^ (4.14) 

we find the correction is more sensitive to the speed of sound Cg. The present ob- 
servation imposes bound on equilateral form of /nl, —256 < /nl < 332, and future 
observation of Planck can detect {J'nlI ^ 5, thus making the correction within the 
sensibility of these observation. 

4.2 K- inflation 

We now consider the correction in K-infiation[20l El] model which also has small 
speed of sound. The Lagrangian of the power law K-infiation is of the form 

P(X,0) = ^i^l(-X + X^), (4.15) 

where 7 is a constant. In the inflationary solution, X remains constant as 



X. = ^. (4.16) 
The solution leads to scale factor a of 

a~t^ (4.17) 
for any < 7 < |. And the speed of sound is 

In order to get small speed of sound, we focus on the region 7 <^ 1. The power 
spectrum in the limit of small 7 is [20] 



1 9 ff2 u 

^ Cs^'yS'K^M:^i ki 

where H is taken to be the Hubble scale at the time of horizon exit for the perturba- 
tions currently at our horizon, and ki is the associated comoving wavenumber. Using 
( I4.18P as well as the fact that data determines ~ 10~^ at horizon crossing, we get 

3 

a72 
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It follows from fl4.20p that the tilt satisfies 

- 1 = -37 + ■ • • , (4.21) 

which allows us to determine 7 ~ ^ and thus ~ ^ using the central value of the 
spectral index in the WMAP results [2]. 

We also find cy = — 1 and cr = in this model, where we choose = 0^. The 
dominant terms in three-point function are those with most Cs in the denominator in 
each kind of coefficients, 

-klk\Kkl){h ■ h)kl{24± + 120^1^ + 720^) 
iH^W^, ^ ^ 1 

-2^-^mnkif - k^,k\kik\){k, . h){8— 

kiko + kok'i + k^ki fcifco/^Sx 
+24 ^5 + 120^) 

+6^) + 2mnkir - kik\k^,k'^,)ki{2^^ 

+Qh±^ + 2A^))+perm.]. (4.22) 
Following the same parametrization as (14.91) . we get 

Ai = -^{{kl)\klf-klk\klkl){k,-k,)^^^^^ 

+ 120^^ + 720^) + perm., (4.23) 

A2 = ^^m)\k',f-kik\k^,e,){h-k,){s—^ 

kiko + kok-i + k-xki kikok-x^ 
+24^^ ^ + 120^^) + perm., (4.24) 

^3 - -^(((^n^+(^?r)(^VA;3)(2i+2Mi±^±^ 
+6^) + 2mmf - kik\kiki)ki{2^^ + Q^-i±^ 

+24 -p^)) + perm.. (4.25) 
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The shape of Ai/kik2k3 is the same as that of DBI case. The shape of A2/kik2k3 
and Azlkik2k^ are drawn in Fig. [2] and Fig. [3] respectively. 




Figure 2: The shape of ^2/fci^2^3- 




Figure 3: The shape of ^3/fcifc2^3- 
In the equilateral triangle configuration ki = k2 = k^ = k, the correction of Jnl 

3 



is divided into IIjlJnlJnl^ 



fl,, = 0.02^, (4.26) 
jl^ = -0.02— (4.27) 

fl^ = 0.25— -f^, (4.28) 
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We estimate the size of (po through Friedmann equation 

SM^i/^ _ ^4 29) 

where 

E = 2XP^x -P = + 3X2), ^4 
and find 0o ~ li^^- Using (14A8D and (ICTj) . we find 

/^^ = -12x/3c^/^^, (4.31) 

and 

/^^ = 5400c^/^^. (4.32) 

Since ~ these two are much smaller than f}^^. Using fl4.18p . fl4.20p and 
= 10^M-\ we find 

fNL ~ 180, (4.33) 

which does not depend on Cs due to the cancelation between and c^. This property 
is different from the case of DBI model. The dominant non-Gaussianity in K-infiation 
without considering noncommutative effect is 

/^^ ^ 0.26c;2 ^ 125. (4.34) 

Since Cs is fixed in K-infiation, it is easy to set an upper bound on Ig. Adding (14.261) 
and (I4.34P together, we get the total non-gaussianity in K-infiation, 

/^i = ISO X IQ-^H^^M^i + 125. (4.35) 

Using the upper bound on the equilateral form of /atl, Inl < 332, we get 

h ^ lO'Mp-,^ (4.36) 



5 Conclusion 

We studied the noncommutative corrections of non-Gaussianities of primordial per- 
turbation in a general framework. The corrections could be large in the models with 
a small speed of sound and a relatively low string scale. We test our result in two 
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particular models, DBI model and K- inflation. We find that the correction of Jnl 
can reach 0(1) or even bigger, and thus is observable within the sensibility of future 
experiments. Our study also shows that the noncommutative corrections are more 
sensitive to the speed of sound than the usual nonGaussianity estimator. This could 
be a clue to distinguish the different contributions to the non-Gaussianity. And also it 
indicates that in the inflation models with large speed of sound, the noncommutative 
correction to the nonGaussianity is small. Moreover, the shape of the corrections are 
different from the commutative case, which can be used as another distinguishing 
estimator. 

From our study, it turns out that the noncommutative corrections become signif- 
icant when the string scale is relatively low. This could be just an illusion, given the 
noncommutative scale may not be string scale. In fact, the leading order correction 
terms are actually proportional to (0^^)^, which could be nontrivially related to string 
scale. If we have a relatively low noncommutative scale, namely we assume a little 
larger noncommutativity, the correction could be larger. In string theory, the non- 
commutative scale depends also on the background B12 field on D-brane. In a sense, 
the corrections tell us the information of the string scale and also of the background 
field. 

On the other hand, without assuming another scale, one can make the natural 
choice (11.51) . Although it is not easy to explain the great difference between the 
Planck scale and the string scale in the perturbative string theory, it is not hard to 
tune the string length scale all the way up to 10^^'^cm in type 1 compactifications and 
nonperturbative heterotic string theory. More interestingly, in recent phenomenology 
study, the string scale is fit by WMAP data to be around 10"^ ~ IQ-'^M^; [811 [821188] . 
Doubtlessly, future experiment will bring forward more data to test this tentative 
choice. 
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A Noncommutative correction of third order per- 
turbation 

The noncommutative correction of third order action of perturbation, which is of 
order 0{e) except the terms proportional to cr, is 



H"^ dl^di^pdi^p + a^di(pdiip + a^d2(pd2(p + a^diipdicp — 2Ha^diipdiip 

-2Ha^d2(pd2(p - 2Ha^diipdiip + {Hafdiipdiip + {Hafd2'^d2'-p 

+ {Haf'diipdiip — dl^pdl^p + 8182^8182^ — {HaafLjp^ — Had(f8f(p 

—Haa(pd2(p + 2Hdl(pd2(p + 2Hd\ipd2ip + 2H^aa(pdf(p + 2H'^aaipd2(p 

—AH 8182^^8182'^ — AH'^8l(p82(p + AH'^8i82(p8i82'p + aa8i^ipCp 

IP 1 

+ad(?2^<^</3 — a^a^ifif — a^(p{28\i~p + 28\i^ — acap)) — - '^" (-^ — l)Haa 

2 00 

(0^ - {8ipf){8lip + 8lip - 2H8fip - 2H8lip + 2Haaip) 



-Pxo'^^^^^^^i^ - imd'i^di^ - ^0,82^8,82^ + (ad)V' 



—aatf^dl^p — aatp^dl^p) — -—^Haa<j)o(p'^{dl(p + d^Lp — 2Hdl(p — 2Hd2ip 



J 5 ^ (pQ Cg zvq Cy a 





(A.l) 



20 



B Noncommutative correction of three-point func- 
tion 



The correction of three-point function of ( is 

[nmk'^r - kimmk, . k,)ki^{2A^^ + 120^^ + 720^) 
+mr + {kin-ih ■ k,)kiki^{2A^^ + 120-^) 

Hk ■ k)kl^{2^, + 6^ + 24^) - 2Aklklklcl^^ 



c 



■s 



+klklc,t{4^ + 12-^) + 12z/c2/c2/c|c,-^ 

+zk,k,c,{2— + Q—) + k^{k2 ■ h)-ij^ + + 12-^) 

r\^^2 ^ ^kik2 + k2k3 + k^ki ^ „kik2k^ 
-ik2 ■ k,)-{- + 2 + 6-^)) 

+{ktk^, + k\kl + k, ■ k,){24tklklklcs^^ - 2tklklc,{2^ + 6-^: 
+kl^c.i^ + ^ + 2^)) - Aiclklklkl± 
+2{{ktmr + {k\f{klf)klkl^{2A^^ + 120-^) 
-Aklk\k^,klklkl^{2A^^ + 120-^) 
+(A;2 • ks)k,ics{- + — ^ + 2-^)] 

+^^?^2^^3^c3(A-12^ + 24^) 
+2zA:,^A:3^c.p?)^ + - - 24^) 

-mlklklcl{±-6^j 

+nc,i/3(i. _ _ klk\klkt){{k, . ^2)^(8-^ 

+24^i^^%±Ml + 120^) - ^?^2^^(24i, + 120^) 
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+mr + {kiy){-{ki ■ k,)ki-{2— + e^— ^ + 24 



, iT r\^i^^ , r,kik2 + k2k3 + k3ki ^ „kik2k3 
. fe)_(2- + 2 + 6-^) 

-(fc^A^S + kikl + h ■ k2)(-klklic,(2-^ + 6||) 
+2*?.e.(l + ^2±^1^2^) 

+2( W)^(fcS)'^ + (t{,2(^.y2)4.1(2_L + ,,h±^ + 24^) 
-4W5*?^(2^+6^5i^ + 24^|^) 

->*Jfc=c,((fc;)= + (k\m± - 6^i^ - 24^) 

•;2;2 3/ ^ /ci + A;2 - 2A;3 kik2 - 2k2k3 - ^kik^ „kik2k3. 
-ik.k^c^i- + 2 — + 6-^) 



2 ' ^4 -/^2L--i-.-.-.^3 

+ {k\f){-?,klklics{2^, + 6-^) 

,r r ^ 1 , nkik2 + k2k3 + kski „kik2k3 . 
. k3-{2- + 2 + 6-^)) 

Off r .ki + k2 kik2. 

-2{k, ■ k2)k,ics{- + — ^ + 2-^) 

+2mmr - kik\kiki)ki^{2^^ + + 24^)] 



-2A;2A;32zc,(2-i^ + 6-^) + iklklkjicl^)] +perm.), (B.l) 

where a and 6 denote the first and second component of k vector respectively, K = 
ki + k2 + ks and perm, denotes all the other terms obtained by rotating the index 
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(1,2,3). All the background value and Hubble constant are calculated at horizon 
crossing. 
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